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Abstract 

We  have  constructed  a  link  between  Physikalisch-Technische  Bundesanstalt  (PTB)  in  Braunschweig,  Germany 
and  the  National  Institute  of  Standards  and  Technology  (NIST)  in  Boulder,  Colorado  using  carrier-phase  Global 
Positioning  System  (GPS)  receivers.  The  link  will  provide  a  direct  measurement  of  the  frequency  difference 
between  UTC(PTB)  and  UTC(NIST),  and  can  also  be  used  to  provide  a  direct  comparison  between  the  primary 
frequency  standards  at  the  two  laboratories.  Based  on  our  previous  work  with  this  method,  we  expect  to  be  able  to 
realize  the  frequency  comparisons  with  an  uncertainty  of  about  2  X  Iff15  using  one  day  of  averaging.  This 
uncertainty  is  smaller  than  the  combined  uncertainties  of  the  primary  frequency  standards  in  both  laboratories, 
and  it  will  therefore  support  a  near-real-time  comparison  of  these  primary  frequency  standards  without  degrading 
their  capabilities  with  the  noise  of  the  transfer  system. 


Introduction 

The  use  of  GPS  carrier-phase  for  time  transfer  on  both  short  (200  m)  and  long  (2400  km)  baselines  has  been 
previously  reported  [1-8].  In  these  experiments  GPS  carrier-phase  time  transfer  has  shown  an  uncertainty  of 
100  ps  and  a  frequency  uncertainty  of  2  parts  in  1015  for  averaging  times  of  a  day.  Initial  short-term 
comparisons  between  two-way  satellite  time  transfer  (TWSTT)  and  GPS  carrier  phase  on  continental 
baselines  have  encouraged  the  use  of  carrier  phase  as  a  feasible  time-transfer  method.  This  paper  reports  on 
an  experiment  to  compare  UTC(PTB)  with  UTC(NIST),  and  CS2  with  NIST7,  the  primary  frequency 
standards  at  Physikalisch-Technische  Bundesanstalt  (PTB)  in  Braunschweig,  Germany  and  the  National 
Institute  of  Standards  and  Technology  (NIST)  in  Boulder,  Colorado,  respectively. 

As  in  previous  experiments,  we  will  estimate  the  relative  clock  behavior  at  6-minute  intervals  using  the  GPS 
carrier-phase  observations.  Details  of  the  geodetic  software  package  and  analysis  are  given  in  [1-4,  9].  As 
independent  estimates  of  the  performance  of  our  analysis,  we  will  use  TWSTT  and  Circular-T  data  between 
NIST  and  PTB.  The  limit  to  the  short-term  performance  of  the  TWSTT  comparison  is  the  nominally  3  times 
per  week  observing  schedule  between  the  two  laboratories.  The  Circular-T  data  are  a  result  of  postprocessed 
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UTC  data  between  various  national  laboratories  and  are  available  monthly,  with  two-three  weeks  latency. 
The  UTC  difference  data  between  laboratories  are  published  at  5-day  intervals  and  primary  frequency  data 
are  published  every  30  days. 

Configuration 

Geodetic  quality  dual-frequency  GPS  receivers  are  installed  at  NIST  and  PTB  (Figure  1).  These  receivers 
produce  both  pseudorange  and  carrier-phase  measurements  at  30  s  intervals. 

The  PTB  GPS  receiver  (PTB1)  is  supplied  with  an  external  5  MHz  reference  signal  from  a  hydrogen  maser 
(H2).  The  primary  standard  at  PTB  is  called  Cesium  2  (CS2).  The  differences  between  UTC(PTB)  and  H2, 
and  CS2  and  H2,  are  reported  on  an  hourly  basis.  H2  is  also  the  reference  for  the  TWSTT  station  at  PTB. 
The  PTB  GPS  receiver  was  installed  in  July  1 999. 

The  NIST  GPS  receiver  (NIS2)  is  supplied  with  an  external  5  MHz  reference  signal  from  an  auxiliary  output 
generator,  AOG2.  The  difference  between  AOG2  and  NIST  hydrogen  maser  one  (Maser  1)  is  reported  eveiy 
two  hours.  The  difference  between  Maser  1  and  NIST7  is  reported  over  a  30-day  interval.  The  output  of 
AOG2  is  steered  to  UTC(NIST),  and  UTC(NIST)  is  the  reference  source  for  TWSTT.  The  NIST  GPS 
receiver  was  installed  in  its  current  location  in  March  1 999;  however,  the  antenna  cable  was  changed  in 
September  1999. 

The  two  receivers  are  approximately  8000  km  apart  and  operated  by  NIST.  The  carrier-phase  data  are 
compared  using  geodetic  software  provided  by  the  Jet  Propulsion  Laboratory  (JPL)  using  other  geodetic 
receivers  that  are  a  part  of  the  International  GPS  Service  (IGS)  network  [9,10]. 

Network  Selection 

In  order  to  make  a  comparison  between  NIST  and  PTB  it  is  important  to  have  other  receivers  in  close 
proximity  to  each  of  the  comparison  sites  to  help  resolve  cycle  ambiguities  (see  Figures  2  and  3).  The 
determination  of  carrier-phase  cycles  is  much  more  difficult  as  the  baseline  length  increases. 

For  this  reason  we  have  selected  a  site  near  NIST  at  Table  Mountain,  Colorado,  TMGO.  TMGO  is  located 
1 5  km  northeast  of  NIST  and  has  a  rubidium  clock  reference. 

We  have  also  performed  the  analysis  with  different  sites  near  PTB:  Potsdam,  Germany  (POTS),  Koetzting, 
Germany  (WTZR),  and  Onsala,  Sweden  (ONSA).  Varying  the  network  sites  allows  us  to  see  the  differences 
in  the  clock  solution.  It  is  important  to  note  that  POTS  and  WTZR  both  operate  receivers  made  by  the  same 
manufacturer  as  those  at  PTB  and  NIST.  ONSA,  however,  operates  a  receiver  from  a  different 
manufacturer.  The  receiver  at  POTS  is  178  km  from  PTB  and  uses  its  internal  clock  as  its  reference.  The 
WTZR  receiver  is  390  km  from  PTB  and  has  a  hydrogen  maser  reference  clock.  ONSA  is  575  km  from  PTB 
and  also  has  a  hydrogen  maser  reference. 
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We  use  intermediate  stations  to  help  define  the  reference  frame  for  the  comparison.  We  have  chosen 
Algonquin,  Canada  (ALGO)  and  the  United  States  Naval  Observatory  (USNO)  in  Washington,  D.C.  ALGO 
is  used  because  of  its  hydrogen  maser  reference  and  well-defined  location  from  Very  Long  Baseline 
Interferometry  (VLBI).  To  help  resolve  cycle  ambiguities  at  ALGO  we  have  chosen  the  USNO  site,  786  km 
away,  which  is  also  referenced  to  a  hydrogen  maser.  The  hydrogen  masers  have  less  noise  in  the  short  term 
and  therefore  provide  better  determination  of  the  correct  carrier-phase  ambiguities. 

Table  1  shows  the  different  network  configurations  that  are  used  for  this  analysis. 

Observed  data  and  Statistics 

Figures  4a-d  show  the  carrier-phase  clock  estimate  solutions  for  the  four  different  networks.  Over  this  17- 
day  time  period  we  compare  carrier-phase  data  with  TWSTT  data.  During  this  time  period  there  are  three 
different  continuous  time  series  of  carrier-phase  data  because  of  jumps  in  the  H2  maser  at  PTB  on  Modified 
Julian  Day  (MJD)  51453  and  MJD  51457.  The  carrier-phase  data  have  been  biased  to  match  the  three 
different  data  series  closest  to  the  TWSTT  data  between  the  jumps.  In  this  plot  we  are  comparing  H2  to 
UTC(NIST).  The  TWSTT  data  are  normally  reported  as  UTC(NIST)-UTC(PTB),  so  a  correction  from 
interpolated  hourly  H2-UTC(PTB)  data  has  been  made  to  the  TWSTT  data  points  in  order  to  make  the 
comparison  with  the  carrier-phase  data.  Over  this  time  period  the  carrier-phase  clock  estimates  appear  to 
agree  with  the  TWSTT  data  in  the  long  term. 

Figure  5  shows  the  fractional  frequency  differences  for  the  UTC(PTB)-UTC(NIST)  carrier-phase  solutions 
from  MJD  5 1 447-5 1 464.  Over  this  1 7-day  time  period  we  compare  the  carrier-phase  data  with  Circular-T 
and  TWSTT  data.  An  hourly  UTC(PTB)-H2  correction  has  been  made  to  the  corresponding  hourly  carrier- 
phase  data  points  to  make  these  comparisons.  For  each  of  the  continuous  time  series  of  carrier-phase  data, 
we  show  the  average  fractional  frequency  across  the  four  network  solutions  as  a  line  covering  the  range  of 
the  data.  The  fractional  frequency  for  each  network  during  that  time  series  is  also  shown.  The  uncertainty 
for  the  5-day  Circular-T  data  is  on  the  order  of  1  part  in  1014.  It  appears  that  the  carrier-phase  solutions 
follow  the  Circular-T  and  T W STT  solutions  with  reasonable  uncertainty,  except  for  near  51457,  where  there 
was  a  reported  jump  in  the  H2  maser  at  PTB.  We  are  investigating  the  source  of  this  error. 

The  fractional  frequency  differences  for  the  CS2-NIST7  comparison  are  plotted  in  Figure  6.  The  Circular-T 
data  for  this  difference  over  this  interval  is  9  parts  in  1 015.  The  first  and  third  series  of  carrier-phase  data 
seem  to  match  the  Circular-T  data  well,  regardless  of  the  network  solution  used;  however,  the  second  series 
of  continuous  data  has  a  frequency  offset  of  approximately  8  parts  in  1015  from  this  solution.  From  our 
initial  analysis  the  uncertainty  on  this  value  is  on  the  order  of  1  part  in  1014,  so  the  frequency  step  might  be 
real. 

Figure  7  shows  the  Allan  Deviation  (ADEV)  for  the  carrier-phase  data  for  MJD  5 1 447-5 1 452.  The  ADEV 
of  the  carrier-phase  data  for  MJDs  51453-51456  and  51458-51464  are  shown  in  Figures  8  and  9, 
respectively.  Only  the  best  and  worst  cases  for  the  four  network  solutions  are  plotted  to  make  the  charts 
easier  to  read.  In  the  time  series  where  we  have  longer  data  spans  the  ADEV  approaches  5  parts  in  1015  at 
one  day. 
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Conclusions 


From  our  initial  analysis  it  appears  that  the  carrier-phase  data  make  it  possible  to  measure  the  differences  of 
CS2-NIST7,  and  UTC(PTB)-UTC(NIST);  however,  there  are  error  sources  that  are  not  yet  fully  understood. 
The  effect  of  the  discontinuities  in  the  carrier-phase  data  combined  with  corrections  made  to  the  data 
requires  further  investigation. 
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Figures  and  tables 


PTB  CONFIGURATION 


TWSTT 


NIST  CONFIGURATION 


Figure  1:  Block  Diagram  of  configurations  at  NIST  and  PTB. 


Figure  2:  Network  Station  Locations  in  the  United  States. 
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Figure  3:  Network  Station  Locations  in  Europe. 


NETWORK  NAME 


STATIONS  IN  NETWORK 


ALGO,NIS2,  PTB1,  TMGO,  USNO.WTZR 


ALGO,NIS2,  POTS,  PTB1,  TMGO,  USNO 


ALGO,NIS2,  PTB1,  ONSA,TMGO,  USNO 


ALGO,NiS2,  PTB1,  TMGO.WTZR 


Table  1:  Network  names  and  descriptions. 


51447  51447.5  51448  51448.5  51449  51449.5  51450  51450.5  51451 

MJD 

--  A  —  B  —  C  —  D  ■  TWSTT  H2-UTC(NIST)  DAoBoCaD 
Figure  4a:  Carrier-phase  and  TWSTT  H2-UTC(NIST)  data.  Symbols  are  added  to  differentiate  the  data. 
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Time  (ns) 


51459.75  51460.25  51460.75  51461.25  51461.75  51462.25  51462.75  51463.25  51463.75 

MJD 

Figure  4d:  Carrier-phase  and  TWSTT  H2-UTC(NIST)  data. 
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Figure  5:  UTC(PTB)-UTC(NIST)  Fractional  Frequency  Difference  data. 
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Figure  6:  CS2-NIST7  Fractional  Frequency  Difference  data. 


Averaging  Time,  x,  (s) 
-*-B  -*-C 


Figure  7:  Allan  Deviation  of  the  carrier-phase  data  for  MJD  51447-51452. 
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Allan  Deviation,  ay(r)  Allan  Deviation,  oy(x) 


Figure  8:  Allan  Deviation  of  the  carrier-phase  data  for  MJD  5 1 453-5 1 457. 


Figure  9:  Allan  Deviation  of  the  carrier-phase  data  for  MJD  5 1 458-5 1 464. 
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Questions  and  Answers 


DEMETRIOS  MATSAKIS  (USNO):  I  just  wondered  if  you  could  comment  more  on  the 
contribution  of  the  clocks  to  the  Allan  deviations  that  you  saw  there.  The  contributions  of  the 
masers,  particularly  the  maser  at  the  PTB,  to  the  Allan  deviations  that  you  showed. 


LISA  NELSON  (NIST):  The  Allan  deviations  that  I  showed  at  the  end,  that’s  the  H-2  to  UTC 
(NIST)  comparison.  You  just  wanted  to  know  what  contribution  the  H-2  makes  to  that? 


MATSAKIS:  Yes,  to  the  clock  itself. 


NELSON:  I’m  unsure  as  to  what  — 


MATSAKIS:  Okay,  we  can  talk  about  it  later. 

DAVID  HOWE  (NIST):  Lisa,  since  you  use  several  stations  to  resolve  the  phase  ambiguity 
issue  and  take  out  those  errors,  can  you  give  us  some  idea  of  how  many  of  those  were  involved 
in  the  whole  data  run? 


NELSON:  Do  you  mean  how  many  biases? 

HOWE:  Right. 

NELSON:  Do  you  mean  the  hardware  resets? 

HOWE:  Yes. 

NELSON:  There  weren’t  any  hardware  resets  other  than  the  fact  that  we  have  the  jump  in 
the  H-2  maser  at  those  two  different  times.  As  far  as  I’m  aware,  the  hardware  was  running 
continuously  through  that  time. 

ED  POWERS  (USNO):  What  about  these  day-boundary  discontinuities?  Do  they  give  you 
much  trouble? 


"K  TT->  T  n  /—V  -k  T 


mars  one  inmg  witn  me  ainerenr  diocks  oi  aaia  mat  we  naa  as  iar  as  aoing  me 
day  boundaries  over  the  6-day  time  period  where  I  actually  had  two  different  series.  Because 
I  do  it  over  a  3-  to  4-day  run  when  I  do  a  run  of  the  different  networks.  We  didn’t  have  any 
problems  as  far  as  that  goes. 


